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Recently, it has been pointed out that two different excesses of events observed at LEP could 
be interpreted as the CP-even Higgs bosons of the MSSM with masses of approximately 98 and 
114 GeV. If this is the case, the entire MSSM Higgs sector is required to be light. In this article, 
we explore such a scenario in detail. We constrain the Higgs and supersymmetric spectrum using B 
physics constraints as well as the magnetic moment of the muon. We then point out the implications 
^ , for neutralino dark matter — next generation direct detection experiments will be sensitive to all 

MSSM models with such a Higgs sector. Finally, we find that all models outside of a very narrow 

| corridor of the parameter space have a charged Higgs boson which will be observed at the LHC. In 

those exceptional models which do not contain an observable charged Higgs, a light top squark will 

' always be seen at the LHC, and likely at the Tevatron. 
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\Q • I. SETTING THE STAGE 

o ; 

The four experiments at LEP have searched for Higgs bosons up to a mass of approximately 115 GeV [1]. Although 
\/np strong indication of a Higgs has been detected by LEP, excesses with statistical significances of 2.3a and 1.7a have 
Q&een reported for Higgs-like events corresponding to masses of 98 GeV and 115 GeV, respectively. Very recently, it 
as been pointed out by Manuel Drees that both of these excesses reported by LEP could be accommodated within 
^e Minimal Supersymmetric Standard Model (MSSM) as signatures of the two CP-even Higgs bosons [2]. 
One of the reason that the LEP excess at 98 GeV has received so little attention might be that the rate observed 
rresponding to this mass is about a factor of ten below that which would be expected from a Standard Model Higgs 
3son. In the MSSM, however, the h-Z-Z coupling is suppressed relative to the value in the Standard Model by a 
lactor of sm(/3 — a), where a is the mixing angle between the two CP-even Higgs bosons and tan/3 is the ratio of the 
# ^vo vacuum expectation values. To predict the number of events seen at LEP for a 98 GeV light Higgs boson, we 
k^mstrain 0.056 < sin 2 (/3 — a) < 0.144 (the precise range of sin 2 (/3 — a), and tuh we allow is the same as used 
^j^L Ref. [2] to allow for easy comparison of results). Meanwhile, the analogous coupling of the heavy CP-even Higgs 
only slightly suppressed below the Standard Model value and is consistent with the number of events reported by 
LEP corresponding to a 115 GeV Higgs boson. These two excesses together yield a combined statistical significance 
of 3. la [2]. 

The mass scale of the MSSM Higgs sector is usually fixed by the mass of the pseudoscalar Higgs boson tua- In the 
well-known leading mf approximation [3, 4], we can very roughly link the pseudoscalar mass to the two scalar Higgs 
masses in the limit of non- mixing top squarks: 

m h + m H m A + m l + 6 
2 2,2 

m H± =m A + m w 

3G F m\ m? 

6 r- o log 7y (1) 

V2ir 2 sin 2 /? m 2 w 

In this approximation is the top squark mass, vn^^A.H are the neutral Higgs masses and ra#± is the charged 
Higgs mass. We immediately see that the existence of two light scalar Higgs bosons leads to all of the MSSM Higgs 
bosons being fairly light. This can lead to interesting phenomenology in the MSSM. In this article, we consider the 
constraints on this scenario from the B — > X s j and B s — > /i + \i~ branching fractions and the magnetic moment of the 
muon. We discuss the implications of such a model on the phenomenology of neutralino dark matter, and find that 
the prospects for the direct detection of neutralinos in elastic scattering experiments to be excellent. We then study 
the spectra of the superpartners and the Higgs sector with two light Higgs scalars and discuss in detail the prospects 
for discovering the charged Higgs boson at the LHC. 
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Figure 1: The mass of the pseudoscalar Higgs versus the charged Higgs is shown in the left frame. In the right frame, /i 
is shown with the charged Higgs mass. For all points shown, 95 GeV < mh < 101 GeV, 111 GeV < ulh < 119 GeV and 
0.056 <> sin 2 (/3 — a) <> 0.144, corresponding to the range matching the observations at LEP. The black points are consistent 
with measurements of B — > X s j at the 3cr level and do not violate the Tevatron constraint on the B s — > branching 
fraction. Blue points violate B s — > but are consistent with B — > X s 7. Green points violate S — > X s 7. These constraints 

will be discussed further in section II. 

In order to study the phenomenology of supersymmetric models with light Higgs scalars, we first performed a scan 
over the relevant parameters of the MSSM. We have varied all masses up to 6 TeV, and tan (3 between 1 and 60. We 
did not assume any specific super symmetry breaking scenario or unification scheme. To take into account the radiative 
corrections to the Higgs masses, we have used FeynHiggs [5], which includes all contributions up to the two- loop level. 
We find that this level of precision is needed to obtain an accurate representation of the MSSM phenomenology within 
this class of models. We give a detailed analysis of the Higgs sector in Section IV. 

In Fig. 1, we confirm the conclusion of Ref. [2] that the entire MSSM Higgs sector is required to be rather light 
to accommodate both excesses reported by LEP. In the overwhelming majority of models found, the mass of the 
pseudoscalar Higgs bosons is between 95 and 130 GeV and the mass of the charged Higgs lies between 110 and 
150 GeV. Although we do find some points outside of this range, in particular a trail of models extending to the upper 
right of the figure (left frame), these are quite rare compared to the dense collection of models near txia ~ 110 GeV 
and m H ± ~ 130 GeV. For each point shown, we require 95 GeV < mh < 101 GeV, 111 GeV < hih < 119 GeV and 
0.056 < sin 2 (/3 — a) < 0.144, corresponding to the range needed to match the observations at LEP. The differences 
between this figure and the corresponding one in Ref. [2] come from the fact that we scan continuously over tan (3. 
In the right frame of Fig. 1, we check the mass of the charged Higgs boson for correlations with the value of \i. We 
have checked that indeed the trails of scattered points at large pseudoscalar masses correspond to the larger values 
of tan/5 the larger the corresponding charged Higgs masses become. Without imposing additional constraints we find 
relatively little correlation, for example, between the Higgs masses and the \i parameter, apart from the chargino mass 
limit of 104 GeV from LEP. 

II. INDIRECT CONSTRAINTS FROM B X s7 , B s fi~ AND (g - 2) M 

Constraints from rare 5-decays and the magnetic moment of the muon are some of the most useful tools we 
currently have to guide our studies of supersymmetric phenomenology. The B physics constraints gain their power 
mostly because of a possible tan (3 enhancement of the bottom Yukawa couplings. In this section, we consider some 
of these constraints within the context of models with 98 and 115 GeV Higgs bosons. 

The branching fraction of B decays to a strange state plus a photon has been measured by the BELLE [6], 
CLEO [7] and ALEPH [8] experiments. The weighted average of these experiments' results indicate BR(£> — > X s j) = 
(3.25 ±0.37) x 10 -4 . In comparison, the Standard Model prediction for this transition rate is (3.70 ±0.30) x 10 -4 [9]. 

The supersymmetric processes which are most likely to contribute substantially to this branching ratio involve a 
charged Higgs or a chargino. These contributions are enhanced by powers of the Yukawa coupling mfrtan/3 for large 
values of tan (3. In passing we emphasize that there are additional A5 corrections [10] which can have huge effects if 
the MSSM spectrum is split between light gluinos or higgsinos and heavier sbottoms, but we will see that this is not 
the part of parameter space in which we are interested in. If the gluino mass is smaller than 300 GeV, the LHC will 
be swamped by gluino pair production with cross sections as large as 1000 pb and the effect of the gluino mass on 
the charged Higgs boson is negligible. 
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Figure 2: The b — > S7 branching fraction as a function of tan /3, the lightest chargino mass, and the charged Higgs mass. Shown 
as horizontal dashed lines are the 2a confidence bounds on this quantity as measured by BELLE, CLEO and ALEPH. Again, 
all points shown have 95 GeV < m h < 101 GeV, 111 GeV < m H < 119 GeV and 0.056 £ sin 2 (/3 - a) £ 0.144. Black points do 
not violate the Tevatron constraint on the B s — > branching fraction. Green points violate this constraint. 

In Fig. 2, we plot BR(£? — > X s 7) versus tan /3 (left frame), the lightest chargino mass (center frame) and the charged 
Higgs mass (right frame). Shown as horizontal dashed lines are the 2a confidence bounds on this branching fraction. 
Again, all points shown have 95 GeV < m h < 101 GeV, 111 GeV < m H < 119 GeV and 0.056 < sin 2 (/3 - a) < 0.144. 
Strictly applying a 2a limit would rule out all points with tan/5 < 10, but on the other hand allowing for a 3a 
deviation from the Standard Model brings all values of tan (3 back into the allowed region. Note that the preference 
of larger values of tan (3 is due to a slight bias toward a finite MSSM contribution: if we force the Higgs sector to be 
light, the 2HDM diagrams will lead to an increase of the observable BR(£? — » X s ^) by typically tens of percent up to 
factor of two. This charged Higgs contributions consist of a tan (3 suppressed term and a constant term, but does not 
exhibit any tan/3 enhancement (if we do not consider anomalously large gluino loops). Because the measured value 
of BR(B — > X s j) is actually slightly smaller than the Standard Model prediction, the chargino has to compensate 
for the 2HDM contribution. Because the chargino contribution to BR(£? — > X s ^) is enhanced by one power of tan (3 
on the amplitude level, we can achieve this by choosing large values of tan/3, as we see in the left panel of Fig. 2. 
We can, of course, try to increase the chargino diagram by making the chargino light, but this is much less efficient, 
because the loop involved is a chargino-stop loop, so that just making the chargino light has comparably little effect. 
Again, we see in the center frame of Fig. 2 that very light chargino masses serve this purpose. Moreover, the large 
chargino contribution has to come with the right sign and therefore prefers positive /i, as is shown in Fig. 1. As we 
can see, the charged Higgs mass and therefore the contribution to BR(£? — > X s j) is basically fixed by the two light 
Higgs masses. Thus merely shifting the charged Higgs mass between 120 and 170 GeV has even less of an effect. 

Limits on the branching fraction of the rare decay B s — > can also be exploited to limit the allowed MSSM 

parameter space. The Standard Model prediction for this branching ratio is BR=(2.4 ± 0.5) x 10 -9 . In particular, 
the Higgs-induced MSSM contributions to this branching fraction scales with tan 6 /3/m^ and tends to be especially 
large in our scenario with a light Higgs sector [11]. If we require the Higgs sector to be light, models with very large 
values of tan/3 are likely to violate the constraints on this quantity placed at the Tevatron [12]. In Fig. 3, we plot 
BR(£? S — > /i + versus tan/3. Shown as a dashed horizontal line is the 90% upper confidence bound placed by the 
Tevatron experiments. The Standard Model value is clearly visible as the low tan (3 nose in the distribution of the 
MSSM parameter points — small tan (3 means small contributions from the MSSM Higgs sector. The upper limit from 
the Tevatron becomes an issue for tan/3 > 10 but it really only develops seriously destructive power for tan/3 > 50. 

Combining the upper limit on the B s — > decay and the measurement of B — > 57 makes it considerably harder 

for us to find viable MSSM scenarios with two light Higgs scalars. On the one hand, we need an MSSM contribution to 
B — » X s 7 from the chargino sector with the right sign and magnitude to compensate for the charged Higgs diagrams. 
This is possible by exploiting the enhancement by one power of tan (3 in the chargino amplitude. On the other hand, 
we do not want too large a contribution to B s — > from the same sector which scales as tan 6 (3. Putting both of 

these constraints together, we favor fairly large, but not too large, values of tan (3 and at the same time accept some 
light superpartners to accommodate B — > X s j. We should stress, however, that this argument is not strict in the 
sense that these two constraints guarantee a light MSSM mass spectrum. By choosing properly tuned parameters, 
we can still get by with a spectrum where the lightest of the charginos and stops weighs more than 1 TeV. 

Looking into the near future, we see that of course the LHC prospects of seeing these light-Higgs MSSM scenarios 
in B physics are excellent. Simulations of CMS, ATLAS and LHCb events probing BR(£? S — > /i + predict tens of 
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Figure 3: Left: The B s — > branching fraction versus tan (3. The dashed horizontal line is the 90% upper confidence bound 

placed on this branching fraction by the Tevatron experiments. Black points are consistent with measurements of B — > X s j 
at the 3cr level. Green points violate this constraint. Right two panels: the contribution to the magnetic moment of the 
muon. The dashed horizontal lines are the 2a bounds from e + e~ data. The black points are consistent with measurements 
of B — > Xsj at the 3a level and do not violate the Tevatron constraint on the B s — >- branching fraction. Blue 

points do violate B s — > but are consistent with £? — > X s 7. Green points violate — > All points shown have 

95 GeV < m h < 101 GeV, 111 GeV < m H < 119 GeV and 0.056 < sin 2 (/3 - a) < 0.144. 

events for 10 fb _1 and the Standard Model decay rate. Moreover, it has been shown that high luminosity triggering 
on this search channel is possible, so that the LHC reach should even cover smaller branching fractions than predicted 
in the Standard Model. Basically all parameter points shown in Fig. 3 will be clearly visible [13]. The only question 
is whether the theoretical and experimental errors will allow us to distinguish between the Standard Model and the 
MSSM predictions. 

Finally, the magnetic moment of the muon has been measured to be anomalously high in comparison to the 
Standard Model prediction. Using e + e~ data, the measured value exceeds the theoretical prediction by 5a^(e + e~) = 
23.9 ± 7.2had-io =t 3.5ibi ± 6 ex p x 10 -10 , where the error bars correspond to theoretical uncertainties in the leading 
order hadronic and the hadronic light-by-light contributions as well as from experimental contributions [14, 16]. This 
measured value is 2.4a above the Standard Model prediction. Experiments using r data, on the other hand, find 
5a^(r + T~) = 7.6 ± 5.8h a d-io ± 3.5ibi ± 6 ex p x 10 -10 , which is only 0.9<r above the Standard Model prediction and 
likely not in agreement with most recent KLOE data [15]. Given this conflict and the marginal statistical significance 
of these measurements, we do not require all of our scenarios to produce the measured value of (g — 2) M . 

Contributions to (g — 2) M occur at the one and two loop levels from diagrams involving both Higgs bosons [17] and 
superpartners [18]. At the one- loop level the contribution from superpartner exchange is proportional to 
tan/3/msusY sign(/x) [18]. The exchange of a light pseudoscalar Higgs contributes like tan 2 (3/m\ in the 

leading order of tan/3. Because we are only considering models with a light Higgs sector, the effects from the A 
exchange will dominate. This is what we see in the tan/3 dependence of the permitted parameters points in Fig. 3. 
In a way, the situation is the same as for B — > Xgj: both measurements sit very slightly away from their respective 
Standard Model predictions (and from a statistical point of view are not very convincing), and the central values can 
be accommodated by choosing large tan (3. In the two right panels of Fig. 3 we see that the B — > X s j constraint 
goes a long way to also accommodate the (g — 2) M measurement in our MSSM parameter space. This is particularly 
striking when we look at the behavior of the points around 5a M = 0. The B — > X s j constraint disfavors the 'wrong' 
sign of Sa^ already, so that the impact of the (g — 2) M measurement on our light-Higgs models is very limited once we 
allow slightly more than a 2a window. Again, we see how very large values of tan/3 are disfavored by B s — > 
which drives the MSSM parameter points toward lighter superpartners, i.e. lighter stops and charginos/neutralinos. 
This constraint is expected to become more stringent over the coming years, even before the LHC will start operation. 

To calculate BR(B —> X s j) : BR(B S — > and 8a ^ we have used the micrOMEGAs program [19]. 

III. IMPLICATIONS FOR NEUTRALINO DARK MATTER 

One very attractive feature of R-parity conserving supersymmetry is that it naturally provides a stable particle 
which, in many models, can be a viable candidate for dark matter. The lightest neutralino is particularly appealing 
in this respect [20] . 
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Figure 4: The thermal relic abundance of neutralino dark matter. The dashed horizontal line represents the density measured 
by WMAP [23]. In the right frame, we compare the relic density to the quantity tan 2 (3 sin 2 6wm 2 z / fi 2 , which for models with 
a bino-like lightest neutralino is a good approximation for e 2 tan 2 /3, and therefore scales with the annihilation cross section to 
bb via pseudoscalar Higgs exchange. The two solid lines represent the analytic result of this calculation using vtla = 120 GeV 
and m x o =50 (bottom) and 500 (top) GeV. All models shown have 95 GeV < m h < 101 GeV, 111 GeV < m H < 119 GeV and 
0.056 < sin 2 (^ — a) < 0.144. The black points are consistent with measurements of B — > X s j at the 3a level and do not 
violate the Tevatron constraint on the B s — ► branching fraction. Blue points violate B s — > but are consistent 

with B — > X s j. Green points violate B — ► X s j. 



Neutralinos can annihilate through a variety of channels, including through the exchange of CP-even or odd Higgs 
bosons, charginos, neutralinos, sfermions and gauge bosons. Which annihilation channel(s) dominates varies from 
model to model. In the scenario we are studying here, however, the presence of light Higgs bosons suggest that 
s-channel Higgs exchange to fermion pairs is likely to be a particularly efficient annihilation channel, especially in 
those models with moderate to large values of tan (3. The formulae for these annihilation channels are collected in the 
appendix. 

The annihilation cross section can be used to calculate the thermal relic abundance present today: 



W - 1^ ^ ? ol v (2) 
Mpi Jgl (a + 3b/x FO ) 



where g+ is the number of relativistic degrees of freedom available at freeze out, a and b are the amplitudes given in 
the appendix, and xfo is the value at freeze out: 

8 Try^FO J 

Over the range of cross sections and masses expected for the lightest neutralino, xpo ~ 20 — 30. To yield the density 
of cold dark matter measured by WMAP [23], we thus require a + 3b/x FO ~ 3 x 10~ 26 cm 3 /s. The main contribution to 
a and b are s-channel CP-odd and CP-even Higgs exchange, respectively. From the amplitudes given in the appendix, 
we find a ~ bx tan 2 /?/3, so we conclude that the 6-term induced through CP-even Higgs plays a subdominant role in 
the relic density calculation to the CP-odd Higgs induced a-term. In the left frame of Fig. 4, we plot the relic density 
of the lightest neutralino versus its mass. We have calculated this quantity using the micrOMEGAs program [19]. 

Over the vast majority of super symmetric parameter space, the composition of the lightest neutralino is dominated 
by its bino component, although a small admixture of wino or higgsino is common. For a bino-like neutralino with 
large or moderate tan (3, the annihilation cross section through pseudoscalar Higgs exchange scales as e 2 tan 2 (3. This 
up- type higgsino component can be approximated in most models by e 2 « e 2 B sin 2 0w sin 2 f3 m 2 z //a 2 . With this in 
mind we plot the neutralino relic density versus the quantity tan 2 f3 I /i 2 in the right frame of Fig. 4. The 

two solid lines shown are the analytic results for tua = 120 GeV and m % o =50 (bottom) and 500 (top) GeV in this 
approximation. The majority of models we find fall within this range, suggesting that their annihilation is in fact 
dominated by and s-channel pseudoscalar Higgs exchange to bb. As expected, no models fall above this range. In 
the models which lie below this range, another annihilation mode or modes (such as t-channel sfermion exchange, 
s-channel Z exchange or s-channel pseudoscalar Higgs exchange to top quark pairs) must contribute substantially, thus 
lowering the relic abundance accordingly. Coannihilations between the lightest neutralino and another superpartner 
may also reduce the relic density. 
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Figure 5: The spin-independent neutralino-nucleon elastic scattering cross section is shown versus the neutralino mass. All 
points shown have 95GeV < m h < 101 GeV, HlGeV < m H < 119 GeV and 0.056 £ sin 2 (/3 - a) £ 0.144. The left frame 
includes models with any relic density, while those models shown in the right frame generate an abundance of dark matter 
within a factor of ten of the observed quantity. The dashed curve is the current limit placed by the CDMS experiment. The 
dotted yellow, solid red, dotted blue and solid magenta curves (from top to bottom) are the approximate projected limits of 
GERDA, CDMS, ZEPLIN-MAX and Super-CDMS (phase III), respectively. The black points are consistent with measurements 
of B — > X s j at the 3a level and do not violate the Tevatron constraint on the B s — ► branching fraction. Blue points 

violate B s — > but are consistent with B — ► X s j. Green points violate B — ► X s j. 



From the right frame of Fig. 4, we can infer that for those models which generate the measured relic abundance 
(those which fall along the horizontal dashed line), most appear to annihilate substantially through pseudoscalar 
Higgs exchange to bb. Models with a bino-like neutralino with a small higgsino admixture and moderate to high value 
of tan (3 are often capable of generating the measured density of dark matter. 

The prospects for the direct detection of particle dark matter are quite encouraging for a neutralino in conjunction 
with a light Higgs sector. In many of the models we study here, the spin-independent neutralino-nucleon elastic 
scattering cross section is dominated by the t-channel exchange of CP-even Higgs bosons. The cross section for this 
process is roughly given by 



SG 2 F m 2 z 



rrip + m x o 



(— €b sin 6^ + £w cos 6w) £ 
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1 2 
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(N\qq\N) 
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where the sum is over quark types and (N\qq\ N) are the matrix elements over the nucleonic state. The quantity cos (3 
should be replaced with sin/5 for up- type quarks. e 2 B , and e 2 H are the bino, wino and higgsino fractions of the 



neutralino. In the sum over quark species, the strange quark contribution dominates with m s (N\ss\N) 
For models with ~ 100 GeV Higgs bosons and a bino-like LSP, we estimate 



0.2 GeV. 
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Of course other elastic scattering channels (squark exchange in particular) can also contribute substantially. In 
Fig. 5 we show the spin-independent neutralino-nucleon elastic scattering cross section for models in this scenario. 
Models shown in the right frame each generate an abundance of neutralino dark matter within one order of magnitude 
of the observed density, while those models shown in the left frame may or may not. The dashed curve is the current 
limit placed by the CDMS II experiment [24]. The dotted yellow, solid red, dotted blue and solid magenta curves 
(from top to bottom) are the approximate projected limits of the GERDA [25], CDMS II, ZEPLIN-MAX [26] and 
Super-CDMS (phase III) experiments, respectively. Edelweiss [27] should also reach a sensitivity similar to that of 
CDMS. To calculate these elastic scattering cross sections, we have used the DarkSusy program [28]. 

From Fig. 5, it is obvious that the prospects for direct detection are excellent for this class of supersymmetric 
models with light Higgs sectors which often dominate the cross section. Most of the favored models are within the 
reach of the CDMS and Edelweiss experiments and those remaining models below this level of sensitivity should be 
within the reach of next generation experiments such as ZEPLIN-MAX or Super-CDMS. 

It is straight forward to understand why models with very small elastic scattering cross sections do not appear in 
Fig. 5. If we limit the magnitude of fi to be less than a few TeV for fine tuning reasons, the higgsino fraction of the 
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lightest neutralino is then bounded from below to be e\ > sin 6w sin /3ra|/(3TeV) 2 - 0.0002. From this, we see 
that (j x n cannot be smaller than 10 -9 to 10 -10 pb without being unnaturally large. 

In addition to these excellent prospects for direct dark matter searches, the characteristics of the lightest neutralino 
in this class of models are fairly ideal for the purposes of indirect detection. Since the annihilation cross section 
during freeze out is dominated by the first term in the expansion, (av) = a-\-bx + 0(x 2 ), the annihilation cross section 
relevant for neutralinos annihilating near the galactic center, throughout the galactic halo, or elsewhere of interest to 
indirect detection, is the maximum value consistent with the measured relic abundance, (av) « 3 x 10 -26 cm 3 /s. This 
makes it rather likely that such models with be within the reach of future indirect detection experiments searching 
for neutralino annihilation products in the form of gamma-rays [29] or anti-matter [30]. The prospects for detecting 
such a neutralino in future cosmic positron experiments are especially promising. 



IV. LIGHT MSSM HIGGS BOSONS AT THE LHC 

Before we begin discussing the prospects of observing a light MSSM Higgs sector at the LHC, we should point 
out that it has been known for a long time that one CP-even Higgs scalar in the MSSM is guaranteed to be seen 
in weak-boson-fusion production with a subsequent decay h,H — > rr [31]. The search strategy is identical to the 
Standard Model search for low Higgs masses and has been extensively studied including detector effects [32]. In the 
special case of three light MSSM Higgs scalars, this discovery channel is known to face a challenge: the rr mass 
resolution will not be sufficient to resolve the two CP-even scalars if their mass difference is less than 0(5) GeV - 
they will appear as one sightly wider Higgs resonance. Moreover, the pseudoscalar Higgs only couples to gauge bosons 
through a dimension-5 operator [33], thus it will not be produced in weak boson fusion. The case of three almost 
mass-degenerate light Higgs scalars has been specifically studied, with an emphasis on distinguishing the three mass 
peaks in the decay to muons [34]. The conclusion is that for Higgs masses below 140 GeV, the inclusive Higgs search 
is challenging and that the heavy Higgs bosons will probably not be separately observable in this channel. The more 
promising strategy is to look for bottom-Higgs associated production with a subsequent Higgs decay to muons and 
at least one b tag. The Yukawa coupling for the heavy scalar and the pseudoscalar Higgs can then be tan (3 enhanced 
and the mass peaks might be resolved for tan (3 > 30 and tua ^ 130 GeV. 

The general lesson we learn from detailed studies by ATLAS and CMS is that it might well be easier to discover a 
charged Higgs boson than the additional neutral scalars. A charged Higgs can be searched for either in anomalous top 
decays, or it can be directly produced as bg — > tH~ . In both cases the decay to a tau lepton is most promising [35]. 
The coverage of models in the (tua - tan (3) plane from this channel has historically included a hole for charged 
Higgs masses just above the top mass. We emphasize that this hole does not mean that the charged Higgs will be 
missed at the LHC in this parameter range. It just means that we will have to combine the search strategies for 
(off-shell) anomalous top decays with the usual associated production [36]. For the purpose of this paper we will use 
a preliminary parton level analysis [37]. We include this result in Fig. 6 

A. Light charged Higgs bosons 

Phenomenologically speaking, the easiest searches for charged Higgs bosons are in the mass range m H ± <m t — vn^. 
A light charged Higgs will almost exclusively decay to a tau lepton, because the decays to Wh° and tb are kinematically 
closed in our scenario. This means we can look for top decays to blv, where there are too many final state taus as 
compared to muons and electrons. In top pair production we let the second top decay into Wb, where the W can 
decay leptonically and hadronically (the latter channel dominates the reach in the charged Higgs mass) [38] . The tau 
from the charged Higgs typically forms a jet. The partial width of a top decay to a charged Higgs contains two terms: 
m 2 tan 2 /3 + m 2 / tan 2 (3. Hence, there will be large branching fractions both for small and large values of tan (3 with the 
weakest point around tan (3 ~ ^mt/rrib ~ 7. Because at the LHC these searches are sensitive to anomalous branching 
fractions of the order of 1% [40], we can expect them to probe charged Higgs masses up to m H ± < m t — mb ~ 160 GeV. 
Including off-shell effects, this reach might even be extended by a few GeV. At the Tevatron the same searches are 
currently in progress, but are strongly statistics limited [39]. They only cover anomalous branching fractions of the 
order of 60%, corresponding to an enhanced Yukawa coupling tan/3 > 20 and masses m H ± < 140 GeV. 

In Fig. 6 we overlay our MSSM parameter points with the LHC (ATLAS) 5a reach. The magenta line corresponds 
to the low-luminosity results from anomalous top decays. The slight bending comes from the mass relation between 
ttia and m H ± as well as from the coupling suppression for intermediate tan/3. We see that most of the parameter 
points with two light Higgs scalars will be clearly visible in anomalous top decays at the LHC, possibly even at the 
Tevatron. 
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Figure 6: Left: The (5<r, 300 fb _1 ) discovery potential for a charged Higgs boson at the LHC. Below the nearly horizontal 
(magenta) line, charged Higgs bosons can be identified through the t — > H ± b, H ± — > rv channel. Models falling to the 
upper right of the red line can be discovered through the channel bg — > tH ± , H ± — > rv. Models with very light A, H ± or 
very small tan/3 are excluded by LEP. Right: the same points shown in the m H + vs tan/3 plane. All models shown have 
95GeV < m h < 101 GeV, HlGeV < m H < 119 GeV and 0.056 < sin 2 (/3 - a) < 0.144. The black points are consistent with 
measurements of B — > X s j at the 3a level and do not violate the Tevatron constraint on the B s — > branching fraction. 

Blue points violate B s — > but are consistent with £ — > X s 7. Green points violate — > X s 7. 



B. Heavy charged Higgs bosons with large tan/3 

The best-studied regime for charged Higgs searches at the LHC is the region of large tan/3 > 10 and charged 
Higgs masses well above the top threshold. Following the approximation in Eq.(l), we see that the required two 
light scalar Higgs bosons imply \/m 2 A + m 2 z + e ~ 150 GeV. Requiring charged Higgs masses above the top threshold 
translates into ttla > 150 GeV, or in other words m 2 z ~ — e. This could be possible if one stop were much lighter than 
the top. On the other hand, the e approximation is based on the assumption of non-mixing stops, so we postpone 
this light-stop scenario to the next subsection. The other way of enhancing the mass difference between the average 
CP-even Higgs mass and the charged Higgs mass (beyond the simple m\ approximation) and to escape the anomalous 
top decays is to create large bottom Yukawa corrections through large values of tan (3 (while keeping the \i parameter 
small). This leads to additional terms contributing to the charged Higgs mass squared, the leading terms of which 
are proportional to raf ra^. The two kinds of contributions which appear with this mixed Yukawa coupling are either 
proportional to log(ra£ ra^ /raf), or they come without this logarithm and are directly proportional to the average 
stop mass instead [4]. We see their effect in Fig. 6 when we compare the tan/3 dependence of tua and ra#±. For all 
value of tan (3 the requirement of two light Higgs scalar limits the pseudoscalar mass to a narrow corridor. We find 
hardly any SUSY scenarios with tua > 110 GeV. By considering large values of tan/3, we can increase the charged 
Higgs mass by 50 GeV and avoid the top decay threshold. 

Because large values of tan (3 are required to at the same time have two light Higgs scalars and a heavier charged 
Higgs, we are automatically driven into the region where the charged Higgs can be found in tH~ production at the 
LHC. The only issue is the hole which usually appears in the LHC coverage region for charged Higgs masses between 
160 and 200 GeV. As mentioned above, this hole is an artifact of the two search strategies meeting in this mass range. 
The associated production process and the anomalous (off-shell) top decay have to be combined to cover this hole. 
This kind of study is ongoing in CMS and in ATLAS [37] and we do not expect this to be a problem once these 
analyses are actually performed on data. In Fig. 6 the densities of points suggests that it is very unlikely to find such 
a light-Higgs MSSM scenario where the LHC sees only one neutral Higgs scalar and no charged Higgs (we find only 
one point somewhat close to this region in our scan). However, we will study these kinds of parameter points and 
their discovery prospects at the LHC in more detail next, to ensure there is no hole in the LHC discovery range. 



C. Heavy charged Higgs bosons and light stops 



As described above, there is this remaining part of the light-scalars parameter space where charged Higgs searches 
in anomalous top decays and in the associated production with a top are challenging. It requires charged Higgs masses 
just above the top threshold and small tan/3. Note that we do not claim that LHC will not see a charged Higgs in 
this region of parameter space. Instead, we emphasize that covering this hole is a crucial task for the near future. 
We perform a dedicated scan around the one point which appears in Fig. 6 close to the hole in the LHC reach for 
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Figure 7: Left: same as in Fig. 6, but with parameters varied as described in the text of section IV C. Right two panels: the light 
stop masses for the same set of parameter points. All models shown have 95 GeV < mh < 101 GeV, 111 GeV < ulh < 119 GeV 
and 0.056 <J sin 2 (/3 — a) <J 0.144. The black points are consistent with measurements of B — > X s ^ at the 3a level and do not 
violate the Tevatron constraint on the B s — > branching fraction. Blue points violate B s — > but are consistent 

with B — > Xsj. Green points violate B — > X s j. 



intermediate tan/3 < 20 and m H ± ~ 180 GeV. We vary each mass parameter up to 50% above or below the value 
of the single point at tan/3 « 8, rriA ~ 140 shown in Fig. 6. Only tan/3 we vary over the entire range 1-60. At this 
specific point, M2 and M3 are both ~ 1 TeV, /i ~ 2 TeV, At ~ —500 GeV, and the bino mass is somewhat light, 
Mi « 140 GeV. 

To see what kind of MSSM parameter choice can enhance the charged Higgs mass compared to the CP-even scalar 
masses, we extend the analytical approximation of Eq.(l) to mixing stops [4]. To avoid the LHC limits, we are now 
interested in parameter points with as small as possible tan /3, so we can neglect the bottom Yukawa coupling. Instead, 
we take into account stop mixing. In two limits the expressions are particularly simple: 

(tan/3 -> 1) 

(tan/3 > 1). (6) 

These results point at another way of increasing the charged Higgs mass while keeping two light Higgs scalar s: we 
choose fairly light stops and increase the trilinear coupling, A t , or the 11 parameter. We have checked explicitly that 
for large values of tan/3 the same can be achieved through increasing Ah. However, keeping both stops light and 
at the same time increasing one way or another the stop mixing parameters, A t — \ij tan/3, means that at least the 
lighter stop and the lighter sbottom will be very light. This is precisely what we observe in Fig. 7, where we see that 
the lightest stop mass is always lighter than approximately 340 GeV. 

Depending on the neutralino and chargino masses there are two possible stop decays we have to consider when we 
search for light stops at the Tevatron and at the LHC. If the stop is heavier than the lightest chargino, it will decay 
mostly into bxt • If this decay channel is kinematically forbidden, as the LEP2 limits suggest for the Run I stop 
searches at the Tevatron, the stop has to decay to qxi- 

The production cross sections for light stops are large: at the Tevatron, a 200 GeV stop will be produced with a 
total cross section of approximately 300 fb, and at the LHC a 300 GeV stop will have a production rate around 10 pb. 
This means that at the LHC, although this kind of MSSM parameter point might be hard to find a charged Higgs 
boson, there will be huge numbers of top squarks establishing supersymmetry. Ongoing Tevatron searches with a 
projected Run II reach of well above 200 GeV (dependent on the neutralino and chargino masses) will already severely 
constrain models with light stops. 

One might worry that in these models with lower tan/3 direct dark matter detection might be more difficult. In 
Fig. 8 we show that this is in fact the case, with spin-independent elastic scattering cross sections as low as 10 -10 pb. 
In each of these models, however, there is a rather light neutralino, m x o < 200 GeV, which makes all of these models 
observable by Super-CDMS and many of them observable by ZEPLIN-MAX, even for <j x n ~ 10 -10 pb (see Fig. 5). 
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Figure 8: The spin-independent neutralino-nucleon elastic scattering cross section is shown versus the CP-odd Higgs mass for 
the range of models described in section IV C of the text. For those models in which charged Higgs searches at the LHC will be 
difficult at the LHC (rriA ^ 140 GeV), this cross section is rather small, <j x n ~ 10 -9 — 10 -10 pb. These models each contain 
a rather light neutralino, however, m x o ~ 70 — 200 GeV, making most of them accessible to ZEPLIN-MAX and all of them 
accessible to Super-CDMS (see Fig. 5). All models shown have 95 GeV < m h < 101 GeV, 111 GeV < m H < 119 GeV and 
0.056 < sin 2 (/3 — a) < 0.144. The black points are consistent with measurements of B — > X s j at the 3a level and do not 
violate the Tevatron constraint on the B s — > branching fraction. Blue points violate B s — > but are consistent 

with B — > X s j. Green points violate B — > X s j. 



V. CONCLUSIONS 

In this article, we describe the phenomenology of models within the MSSM containing CP-even Higgs bosons with 
masses of approximately 98 and 114 GeV, motivated by the excesses observed by LEP [2]. We study in detail indirect 
constraints on such light-Higgs models, coming from B — > X s 7, (g — 2) M and B s — > jiji. The first two constraints 
require fairly large tan (3 when compared with data, the latter tends to disfavor this regime, and the upcoming Tevatron 
Run II results are going to close in on our light-Higgs parameter space. However, it is possible to avoid all indirect 
constraints by tuning the different weak-scale SUSY breaking parameters. 

We find that the kind models we are interested in typically include a light charged Higgs boson along with the other 
light Higgs bosons. Such a state is observable in anomalous top decays at the Tevatron or at the LHC. Scenarios with 
charged Higgs above the top threshold are automatically driven in the larger tan (3 regime, where the charged Higgs 
can be detected in bottom-gluon fusion at the LHC. We find exceptionally models which fall between these two mass 
regions, so that neither of the two searches is optimized and a combination of the two search tactics would be needed. 
However, these challenging scenarios require one very light top squark (m^ < 300 GeV), so a large fraction of them 
will already be ruled out at the Tevatron. 

We also explore the phenomenology for neutralino dark matter in this class of models. The prospects for direct 
detection are excellent, with the majority of models within this scenario being testable in currently operating exper- 
iments, such as CDMS-II. We find that all models within this class will be testable in next generation direct dark 
matter detection experiments, such as ZEPLIN-MAX or Super-CDMS. The prospects for indirect are also favorable 
in this scenario. 
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help with FeynHiggs. DH is supported by the Lever hulme Trust. 
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VI. APPENDIX: NEUTRALINO ANNIHILATION CROSS SECTIONS 



The squared amplitudes for the processes, x°X° 
angle are given by [21] 



ri2 r<2 



A — > // and X X ~~ ^ H — > //, averaged over the final state 



Ami 



ff ~ (s - m 2 A ) 2 + m 2 A Y\ 16tt 



J ff 



C ffH C l° x °H ( s ~ ^ m x o)(s - Am 2 ) 



(s-m 2 H ) 2 + m 2 H T 2 H 



16tt 




(7) 



(8) 



Here the label H denotes either CP-even Higgs boson. CffA, C x o x o A , CffH and C x o x o H are the fermion-fermion-Higgs 
and neutralino-neutralino-Higgs couplings. Ta,h are the widths of the respective Higgs bosons. 

These squared amplitudes can be used to attain the thermally averaged annihilation cross section [22] 



(av) 



m x° 



mf 



3T 



;(s) 



2w'(s ) 



T + 0(T 2 ) 
+ 0(T 2 ), 



(9) 



s^4m 2 +6m T 



where T is the temperature. Keeping terms to zeroth and first order in T is sufficient for the relic abundance 
calculation. Writing this as an expansion in x = T/m % o, (av) = a + bx + G(x 2 ), we arrive at 
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\c f m 2 cos 2 a (™' " m /V 1_m // m * 



Wirm^ sin 2 /? (4m 2 -m 2 h + im h T h ) 2 
x [ec/ (ejy — tan^) sin a — cd (tw — £b ta,n$w) cos a] 2 



(10) 



where Cf is a color factor, equal to 3 for quarks and 1 for leptons. The quantities cb,w,u,d are the components of the 
lightest neutralino which are bino, wino, up- type higgsino and down- type higgsino, respectively. These expressions 
assume that the final state fermions are down- type. For up- type fermions the factor of tan 2 /? should be replaced by 
cot 2 ft. 
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